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Abstract 

We analyze r lepton decay observables, namely moments of the hadronic spectral 
density in the finite energy interval (0, M T ), within finite order perturbation theory 
including a A s corrections. The start of the asymptotic growth of the perturbation 
theory series is found at this order in a scheme invariant manner. We establish 
the ultimate accuracy of finite order perturbation theory predictions and discuss the 
construction of optimal observables. 



1 Introduction 



The study of r lepton decays provides a wealth of information on low energy hadronic 
physics where the accuracy of experimental data is permanently improving [I], [| . The 
central quantity of interest is the hadronic spectral density. The spectral density has 
been calculated with a very high degree of accuracy within perturbation theory (e.g. 
||, |], The structure of observables - related to the two-point correlator of hadronic 
currents with well established and simple analytic properties - makes the comparison 
of experimental data with theoretical calculations very clean. All these features make 
r lepton physics an important area of particle phenomenology where theory (QCD) 
can be confronted with experiment to a very high precision || [7[ fs|, P, PH|. 

In the present note we show that within the finite order perturbation theory 
analysis the ultimate theoretical precision has been reached already now. The limit 
of precision exists due to the asymptotic nature of the perturbation theory series. 
The actual magnitude of this limiting precision depends on the numerical value of 
the coupling constant which is the expansion parameter. We perform our analysis 
and reach our conclusions in a renormalization scheme invariant way. 

The normalized r lepton decay rate into hadrons h is given by 

R, = ^^ = N c (l + 6) (1) 

1 [T — > IVV) 

with 

R e T xp = 3.649 ± 0.014 and 5 exp = 0.216 ± 0.005 . (2) 

The first term in eq. ([I]) is the parton model result while the second term 5 represents 
the effects of QCD interaction. In this paper we neglect electroweak corrections 
altogether. The theoretical expression for the rate is given by 

r M T s <? fin 

R T = N C 2(1- ^-) 2 (l + 2^-)p(s)^- . (3) 
c Jo K Ml' y Ml wy 'Ml y ' 
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The spectral density p(s) is related to Adler's D-function through the dispersion 
relation 

The D-function is computable in perturbation theory. In the MS scheme the pertur- 
bation theory expression for the D-function is given by 



D(Q 2 ) = 1 + ^ + 1.64 ( - ) 2 + 6.37 f^V + k 3 ( + . . . (5) 



7T V 7T / V 7T / \ 7T 



where the running coupling is normalized at the scale fi = Q. The light quarks w, <i 
and s are taken to be massless. Eqs. (|3|-|5]) constitute the full theoretical information 
necessary for the perturbation theory analysis of the r system. The fourth order 
MS-scheme coefficient is not known at present. 

In the present note we do not systematically discuss non-perturbative effects stem- 
ming from standard power corrections ||11|| . Also, the infinite resummation of the per- 



turbation theory series different from the standard renormalization group improve- 



ment is used only as a toy example [T3]. The standard power corrections due to 



nonvanishing vacuum expectation values of local operators within operator product 
expansion are relatively small and can be simply accounted for if desired. The coeffi- 
cient functions of local operators are known in low orders of the perturbation theory 
expansions and there is no necessity to thoroughly analyze their convergence prop- 
erties. It is the high precision achieved in the experimental analysis of r decays and 
the rather advanced stage of theoretical description that calls for a detailed analysis 
of the physics of the r system. 
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2 Internal perturbation theory description of basic 
r system observables 

The central quantity of interest in the r system is the hadronic spectral density 
which can be measured in the finite energy interval (0, M T = 1.777 GeV). Being a 
distribution (in theory) or a rapidly varying function in the vicinity of resonances 
(in experiment) the hadronic spectral density cannot be analyzed pointwise within 
perturbation theory. The appropriate quantities to be analyzed are the moments (or 
generalized Fourier components over a chosen complete set of test functions). We 
define moments of the spectral density by (with M T chosen to be the unity of mass) 

M n = (n + 1) f 1 p(s)s n ds = l + m n . (6) 
Jo 

Due to the completeness of the basis {s n : n = 0, . . . , oo} the moments m n contain 
the entire information about p(s). The invariant content of the investigation of the 
spectrum, i.e. independent of any definition of the charge, is the simultaneous analysis 
of all the moments. Note that within finite order perturbation theory the moments 
eq. (^) coincide with the results of contour integration ]TB|, [TJ], [T^, [RJ because of 
analytic properties of the functions ln p s. 

In order to get rid of artificial scheme-dependent constants in the perturbation 
theory expressions for the moments we define an effective coupling a(s) directly on 
the physical cut through the relation 

p(s) = l + a{s). (7) 

All the constants that may appear due to a particular choice of the renormalization 
scheme are absorbed into the definition of the effective charge e.g. fTT| , r8[ |I9 , pO 



Note that if there was no running (as in the conformal limit of QCD with vanishing 
/^-function or at the infrared fixed point) then the whole physics of the r system 



in the massless approximation (without strange particles, for instance, and including 
only perturbative corrections without possible power corrections) would reduce to the 
determination of a single number a(M T ) = a and consequently there would not be 
any problems with the convergence of the perturbation theory series. Because of the 
running of a(s), however, different observables, i.e. different moments of the spectral 
density, generate different perturbation theory series from the original object p(s) in 
eq. (|7|). The whole set of moments needs to be analyzed in a scheme invariant way e.g. 
2T| , P^ , p5| . Note that the introduction of a natural internal coupling parameter such 



as the effective charge a(s), allows one to extend the perturbation theory series needed 
for the description of relations between observables by one more term as compared 



to the analysis in e.g. the MS scheme (e.g. [pi , [24|j). When defining the effective 
charge directly through p(s) itself we get theoretical perturbative corrections to the 
moments only because of running. Without running one would have 

M n = 1 + a(M T ) = 1 + a or m„ = a (8) 

and the perturbation theory analysis would be over (we neglect power corrections for 
the moment!). In any given order of perturbation theory the running of the coupling 
a(s) defined in eq. (0) contains only logarithms of s with coefficients given by an 
effective /3-function 

a(s) = a + p La 2 + (fcL + /3 2 L V + {faL + ^PiPqL 2 + /3 3 L 3 )a 4 + . . . (9) 

where a = a (M T 2 ), L = \n(M 2 /s). The contributions of powers of logarithms to the 
normalized moments are 

(n + l)j\^(l/s)ds = j^- p . (10) 

Therefore, at fixed order of perturbation theory the effects of running die out for large 
n improving the convergence of perturbation theory series. With the definition of the 



charge according to eq. (|7|) all high order corrections vanish at n — > oo at any fixed 
order of perturbation theory. With running one has instead of eq. (^|) 

m = a + 2.25a 2 + 14.13a 3 + 87.66a 4 + (433.3 + 4.5fc 3 )a 5 

m x = a + 1.125a 2 + 4.531a 3 + 6.949a 4 + (-175.2 + 2.25A; 3 )a 5 

m 2 = a + 0.75a 2 + 2.458a 3 - 1.032a 4 + (-142.6 + 1.5k 3 )a 5 

m 3 = a + 0.563a 2 + 1.633a 3 - 2.542a 4 + (-110.4 + 1.125A;3)a 5 

m 100 = a + 0.022a 2 + 0.041a 3 - 0.25a 4 + (-4.08 + 0.045A;3)a 5 (11) 

For large n the moments behave better because the infra-red region of integration is 
suppressed. Note that the coefficients of the series in eq. fllT]) are saturated with the 
lowest power of logarithm for large n for a given order of perturbation theory, i.e. 
they are saturated with the highest coefficient of the effective /3-function. 

Higher moments are not welcome from the experimental point of view. They are 
dominated by the contributions coming from the high energy end of the r decay spec- 
trum (therefore they converge better perturbatively) but experimental accuracy for 
the moments basically deteriorates with increasing n because poorly known contri- 
butions close to the right end of the interval are enhanced. To suppress experimental 
errors from the high energy end of the spectrum the modified system of mixed mo- 
ments 

Mm = £ P (s)(l - s) k s l ds = 1 + m« (12) 

can be used e.g. The weight function (1 — s) h s l has its maximum at s = l/(l + k). 
The integral in eq. ([i~2[ ) is dominated by contributions from around this value. A 
disadvantage of choosing such moments is that the (1 — s) k factor enhances the infra- 
red region strongly and ruins the perturbation theory convergence. As an example 
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one has 



m o = a + 2.25a 2 + 14.13a 3 + 87.66a 4 + (433.3 + 4.5fc 3 )a 5 

m 10 = a + 3.375a 2 + 23.72a 3 + 168.4a 4 + (1042. + 6.75fc 3 )a 5 

m 20 = a + 4.125a 2 + 31.24a 3 + 241.1a 4 + (1683. + 8.25A; 3 )a 5 

m 3Q = a + 4.688a 2 + 37.51a 3 + 307.3a 4 + (2324. + 9.375A;3)a 5 . (13) 



The values of the coefficients in the series eq. ([131) can be found in a concise form 
for arbitrary k at any giving finite order of perturbation theory. For instance, the 
contribution of the log-term is given by 

fc+i 



(Jfe + l) s) k \n(l/s)ds = Y- 

Jo j=l J 



(14) 



In contrast to eq. (|10"D it increases as In (A;) for large k making the coefficients of the 
perturbation theory series large. The contribution of log 2 -term reads 

i /jfe+i i \ 2 fc+i -i 

(k + 1) / (1 - s) k \n 2 (l/s)ds = E - +E- (15) 
J0 \j=i 3 ) j=i 3 

and can be seen to grow as m 2 {k) for large k. 

In practical applications our formal criterion of the accuracy which the series 
provides is given by the numerical magnitude of the last term of the series. However, 
this criterion should be applied with great caution. Because of the freedom of the 
redefinition of the expansion parameter the last term of the series can always be made 
arbitrary small for any given observable. One can give an invariant meaning to the 
quality of the perturbation theory expansion only for a set of observables. 

Before proceeding we would like to comment on the contribution of power correc- 
tions to the systems of moments eqs. (|6|, |l2"l) and the interplay between the magnitude 



of this contribution and the structure of the perturbation theory series. For the sys- 
tem of moments in eq. (|]) the contribution of power corrections reduces to a single 
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term (neglecting the weak log(Q) dependence of coefficient functions of the oper- 
ator product expansion which is a common practice) of the form (A 2 /M^) n which 
decreases very fast with n (A is a typical scale of power corrections related to the 
non-perturbative scale of QCD and A < M T ). This makes the perturbation theory 
contribution dominant in the total result. This perturbative term for the large n 
moments is saturated with high energy contributions and therefore converges pertur- 
batively. The convergence becomes even better with increasing n. The moments are 
perturbatively dominated and, therefore, precise. On the contrary, for the system 



of mixed moments in eq. (12) with I ~ the large k moments are saturated with 
low energy contributions, i.e. basically with the contribution of the ground state res- 
onance, and therefore are completely nonperturbative which is reflected in the fast 
deterioration of perturbative convergence. The contribution of power corrections to 
the moments eq. (0) picks out many terms in all orders from n = 2ton = k + l. 
Nothing definite can be said about such a sum of power corrections in any realistic 
case. This also indicates the importance of power corrections for mixed moments. The 
perturbation theory series is the same both for the vector and axial vector channels 
while the lowest resonance contributions are completely different (the pion instead of 
the p- meson). Therefore no method of summation of perturbation theory series can 
bring it to agree with experiment. In this case perturbation theory is in trouble and 
the power corrections provide the correct result for large k. Large k mixed moments, 
therefore, are not usable within the perturbation theory framework even if they are 
preferable from the experimental point of view. 

Thus, one faces the usual clash between experimental and theoretical accuracy 
which is reflected in our case in the range of (k, I) values for the mixed moments that 
are chosen as optimal observables. Having explicit perturbation theory formulas at 
hand one can establish the ultimate theoretical accuracy implied by the asymptotic 
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character of perturbation theory series for a given experimental observable with any 
stated precision. This allows one to conclude which error - experimental or theoretical 
- dominates the uncertainty of an observable related to r decay physics. 

For our numerical estimates we take a = 0.111 as obtained from the corresponding 
value of the MS charge. From the set of moments {m n ; n — 0, . . . , } the moment niQ 
has the biggest infra-red contribution. Therefore a set of observables has the worst 
accuracy if the moment mo is included in the set. For mo one obtains 

m = 0.111 + 0.0277 + 0.0193 + 0.0133 + (0.0073 + 0.000076fc 3 ) . (16) 

As mentioned before the numerical value of the coefficient fc 3 is unknown at present. 
In some of the following evaluations we want to fix its value to have a feeling of the 
importance of the last term of the perturbation theory series. One popular value is 
&3 = 25 based on Pade approximation. Another value fa A — 91 nullifies the fourth 
order coefficient of the effective f3- function Both these numbers are used only 



for illustrative purposes while our conclusion about perturbatively commensurate 
observables is independent of fc 3 . Numerically for k 3 = 25 one has for the zeroth 
moment 

m = 0.111 + 0.0277 + 0.0193 + 0.0133 + "0.009" . (17) 

Formally, the convergence still persists even in eq. (|T7|) if one only requires subsequent 
terms of the series to decrease but the convergence is very slow. Also, the total 
contribution of the four higher order terms is more than 60% of the leading one. For 
the first moment the convergence is considerably better 

mi = 0.111 + 0.014 + 0.006 + 0.001 + (-0.003 + 0.00004fc 3 ) (18) 

and for fc 3 = 25 

mi = 0.111 + 0.014 + 0.006 + 0.001 - "0.002" . (19) 



With the choice fc 3 = 25 the 0(a 5 ) term already shows numerical growth. If one 
keeps only the smallest term one gets a formal accuracy of about 1% and the total 
contribution of the three higher order terms gives only about 20% of the leading one. 
The large difference in accuracy between m and mi is a general feature of the moment 
observables at fifth order of perturbation theory: one cannot get a uniform smallness 
at this order for several moments at the same time adjusting only one number k%. For 
a = 0.111 we therefore conclude that one is starting to see the onset of asymptotic 
growth at fifth order. The growth of the terms is independent of any definition of 
the charge if several moments are analyzed simultaneously and this feature cannot 
be changed by any choice of k^. For any single moment, e.g. m , one can always 
redefine the charge and make the series converge well at any desired rate but then 
other moments become bad in terms of this charge. The invariant statement about 
the asymptotic growth is that the system of moments m n with n = included cannot 
be treated perturbatively at the fifth order of perturbation theory for the numerical 
value of the expansion parameter a = 0.111 if one wants to obtain an accuracy 
better than 5% - 10%. This statement about the ultimate accuracy of the set of 
moment observables attainable in fifth order of perturbation theory is independent 
of whichever numerical value k 3 takes. If, however, the system of moments m n does 
not include m as an observable, a uniform accuracy better than 1% can be obtained 
for such a system within perturbation theory. For instance, excluding mo and using 
k 3 ~ 100 one can make the system of moments with n > 1 perturbation theory 
commensurate at fifth order in a sense that all fifth order terms can be made small 
simultaneously. To demonstrate this in a scheme invariant way we choose the second 
moment (which is already well convergent) as a definition of our experimental charge 
and find 

m = m 2 + 1.5ml + 9.417m^ + 59.28m2 + (310.3 + 3A; 3 )m^ 
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mi = m 2 + 0.375m2 + 1.51ml + 2-527mj + (-54.45 + 0.75k 3 )ml 

m 2 = m 2 

m 3 = m2 - 0.19ml - 0.544m^ + 0.1A2m\ + (35.2 - 0.375Jfe 3 )77iii 

m 4 = m 2 - 0.3ml ~ 0.803m^ + 1.69m| + (56.641 - 0.6Jfe 3 )mij . (20) 

The convergence for the moments mi-m^ (and for n > 4) is fine. The total contri- 
bution of higher order corrections is small. The worst series is the one for the zeroth 
order moment. Eq. (|20"D shows that no choice of k 3 yields accuracy for both mo and 
mi which is essentially better than the fourth order term. In fact, there is a narrow 
window 40 < k$ < 60 where the formal criterion of convergence is satisfied for both 
m and mi but we do not find it natural to rely on such a fine tuning and even 
then the accuracy of the zero moment is only about 10%. This is an indication that 
the ultimate accuracy of the perturbation theory expansion for the zeroth moment 
has been reached. If the moment mo is excluded the choice k 3 ~ 100 allows one to 
make the convergence fast even to fifth order and no conclusion about an asymptotic 
growth is possible. 

The perturbation theory expansions for the system of moments with (1 — s) h 
weight shows worse behavior. With the above criterion of accuracy, the precision 
which is given by the series from eq. ( |T3] ) is of order 10% - 20% for the numerical 
value of a. This is not enough for a comparison with experiment at the present level 
of precision. For instance, an expansion of the higher moments in eq. ( pT3| ) in terms 
of the first one (which is the most perturbative one for this system) goes as follows 

moo = m o = 0.17 

m 10 = 0.17 + 0.033 + 0.022 + 0.011 + (-0.005 + 0.00032^) 

m 20 = 0.17 + 0.054 + 0.043 + 0.027+ (0.0015 + 0.00053^) 

m 30 = 0.17 + 0.070 + 0.061 + 0.046 + (0.014 + 0.0007^) (21) 
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These series possess a formal accuracy of from 6% to 25% and the contribution of 
higher order terms can be as large as the leading term. Because of the slow conver- 
gence there is no sign of improvement with higher order of the perturbation theory: 
the series expansions do not allow any reliable estimate of accuracy for large mixed 
moments. Also while for the moments eq. (§) the total contribution of corrections 
is small, the situation is different here. The total change of the leading order result 
due to higher order corrections is considerable and strongly differs for various mo- 
ments. This is another indication that the set of mixed moments is not commensurate 
perturbatively. 



3 r decay rate 

The t decay width is given by a specific linear combination of moments. Because 
of the factor (1 — s) 2 present in eq. ([3]) the convergence property of the total decay 
rate observable is not optimal. The (1 — s) 2 factor enhances the infra-red region of 
integration, i.e. the relative magnitude of the contributions of logarithms \n(M 2 /s) 
at small energy. The concrete shape of the weight function with the weight factor 
(1 — s) 2 is the main source of slow convergence. One has 

r T = a + 3.563a 2 + 24.97a 3 + 174.8a 4 + (1041. + 7.125fc 3 )a 5 . (22) 

Using a = 0.111 and k% = 25 

r T = 0.111 + 0.044 + 0.034 + 0.027 + "0.021" (23) 

Formally the consecutive terms decrease but the decrease is very slow. One can see 
that the pattern of convergence mainly follows that of the moment rfi2o from eq. ( |1"3"D 
because of the factor (1 — s) 2 in eq. (|3|). 



12 



Eqs. (|2|) and ( p3| ) show the essence of the problem we are addressing. In the finite 
order perturbation theory analysis one has to compare 6 exp with 6 th 

0.216 ± 0.005 = 5 exp vs 5 th = 0.111 + 0.044 + 0.034 + 0.027 + "0.021" (24) 

and the uncertainty of the theoretical expression "0.021" (or even 0.027) is much 
larger than the experimental error 0.005. Thus, the theoretical uncertainty due to 
the truncation of the perturbation theory series is much larger than the experimental 
error of the corresponding observable. The common practical expectation for theo- 
retical perturbation theory expansions to be useful is the smallness of the total higher 
order corrections if nothing is known about the convergence of the expansion. For 
the rate observable the corrections increase the leading order result by a factor of 2. 
Note that one can improve the explicit convergence of the rate observable by a spe- 
cial redefinition of the expansion parameter due to renormalization scheme freedom. 
However, then the first moment of the differential decay rate will behave wildly. It 
is this feature that prompts us to reach definite conclusions about the asymptotic 
growth of perturbation theory expansion independent of any scheme. Two different 
sets of observables where one set includes the moment m and the other set does not 
include it are not perturbatively connected with an accuracy required by experiment. 
Indeed, the first s-moment of the differential decay rate dR T /ds gives the series with 



faster convergence than eq. fl22 ) 



rW = a + 2.138a 2 + 10.15a 3 + 28.43a 4 + (-268.3 + 4.275A; 3 )a 5 (25) 

or numerically with /c 3 = 25 

rf } = 0.111 + 0.026 + 0.014 + 0.004 - "0.003" . (26) 

The second s-moment has even better perturbative expansion 

r ( 2 ) = a + 1.575a 2 + 6.186a 3 + 6.386a 4 + (-283.3 + 3.15£;3)a 5 (27) 
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and numerically with k% = 25 

rf } = 0.111 + 0.0194 + 0.0085 + 0.001 - "0.003" . (28) 

The fifth order term is larger than the fourth order term for k 3 = 25. No choice of 
k 3 can simultaneously make all these three observables convergent at fifth order. If 
one chooses k% ~ 100 in order to guarantee for a better convergence of the higher 
moments (which is physically motivated) one almost destroys the perturbation theory 
series for the decay rate eq. (|22|) . 

The (1 — s) n moments of the differential decay rate suppress poorly known high 
energy experimental data. Taking n = 1 as an example one has 

rf"*) = a + 4.173a 2 + 31.31a 3 + 237.6a 4 + (1603. + 8.35/c 3 )a 5 (29) 

and numerically for k 3 = 25 

rf ~ s) = 0.111 + 0.051 + 0.043 + 0.036 + "0.031" . (30) 

For &3 = 100 the series reads 

rf ~ s) = 0.111 + 0.051 + 0.043 + 0.036 + "0.041" (31) 

which gives only about 30% accuracy and more than a factor 2 change of the leading 
order term. We conclude that the theoretical precision cannot compete with the 
experimental precision. 

There are two distinct problems in analyzing r decays: one is to describe the 
set of observables of the system using its internal coupling parameter defined to get 
the highest precision and establish whether the set is perturbatively commensurate, 
while another is to extract the standard MS parameters. It can happen that the set of 
observables is perturbatively connected with some given accuracy but the MS coupling 
a s is not the best parameter for the expansion. This is the case here. In internal terms 
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the r system is described with higher accuracy in terms of the number of perturbation 
theory terms than in the MS scheme. However, at this level of expansion one sees the 
asymptotic growth of the perturbation series for the numerical value of the expansion 
parameter fixed by experiment. 

The expression for the decay rate in the MS scheme possesses only O(af) accuracy 

» = (v) + 5 - 20 (v) 2 + 2fM (7)' + (78 -° + h) (v) 4 (32) 

with a numerical precision of only 30% again. A numerical value for a s is usually 



extracted treating the three first terms of the rate expression eq. (32) as an exact 
function. The numerical value found is rather precise. However, the accuracy of 
the numerical prediction for other observables is dominated by the uncertainty of 
truncation of the series and is poor if the observable contains the zeroth order moment. 
Therefore the comparison of different observables of the system cannot be done with 
high precision and the ultimate precision is limited by the asymtotic growth of the 
perturbation theory series. The coupling constant, though important, is still an 
artificial parameter and the knowledge of its precise numerical value does not suffice 
for computing observables with sufficiently high precision. 

The investigation of the r system can be performed in N 4 LO without any free 
parameters with the use of the internal charge a (and even N 5 LO with the single free 
parameter k% which does not affect the conclusion about the asymptotic structure 
of perturbation theory series). However the MS scheme coupling can be expressed 
through a only up to NNLO because of the unknown coefficient k 3 . The extraction 
of the MS charge from a can be made through the relation 
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a - 1.64a 2 + 15.7a d + (49.6 - k 3 )a 4 + ... (33) 



with the reasonably fast convergence for k 3 = 25 or k 3 = 100. 
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4 Infra-red fixed point as a model for an infinite 
perturbation theory series 

The accuracy of approximating a function with the sum of a finite number of term 
of its asymptotic series depends strongly on the analytic structure of this function. 
Generally, there is an infinite number of ways to sum an asymptotic series with 
quite different results. Therefore estimates of accuracy based on asymptotic series 
along can be rather misleading. To discuss this issue in more detail (though for 
illustrative purposes only) we consider a model for the exact function as source for 
the perturbative expansion (or recipe of the infinite resummation) . The model uses 
the existence of the infra-red fixed point for the running coupling in the third order 
of perturbation theory which allows one to extrapolate running to the origin. In this 
particular case we can compare perturbative expansions with an exact answer. This 
example allows one to check the general conclusions about the asymptotic structure 
and the divergence of the series even if this is done in a model dependent way. The 
effective /3-function is given by 

P eff (a) = -^a 2 - 4a 3 + 25.7a 4 + (409.5 - ^k 3 )a 5 + 0(a 6 ) (34) 

where the only free parameter is k 3 because the (3% coefficient in the MS scheme is 
known ||26|| . The third order approximation of the /3-function eq. (|34|) possesses an 
infra-red fixed point with the value cif = a(0) = 0.384. 

The integration in eq. @ can be explicitly performed with the third order f3- 
function from eq. (j34"t). With the initial value a = 0.111 one obtains ml = 0.1605, 
m{ = 0.130954 to be compared with the results of eqs. (p| ,l^,p!8). The naive estimate 
of the accuracy does not always work for this resummation recipe for all the moments. 
Indeed, having the explicit model at hand one can generate an arbitrary number of 
terms of the perturbative expansion. For the zeroth moment the series diverges with 
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the pattern 

m{ ix = 0.111 + 0.028 + 0.019 + 0.013 + 0.014 + 0.018 + 0.029 + 0.053 + 0.114 + .. . (35) 

giving an ultimate accuracy of only about 10%. This accuracy is obtained keeping 
the smallest term. The sum of the first four terms gives the best estimate 

m fix,best = qui + Q.028 + 0.019 + 0.013 = 0.171 ± 0.013 

to be compared with the exact result ml = 0.1605. The central value is a bit too 
high but still within the uncertainty given by the last term. 

For m[ tx one also finds a divergent series but with a much faster decrease of the 
first few terms. The pattern of "convergence" is given by the following huge expression 
which we display to show how complicated things can become. One has 

m{ ix = 0.111 + 0.013861 + 0.006197 + 0.001054 + 0.000480 + 0.000088 

+ 0.000053 + 0.000016 + 0.000015 + 0.000014 + 0.000019 + 0.000026 

+ 0.000042 + 0.000072 + 0.000135 + 0.000268 + 0.000568 + 0.001277 

+ 0.003 + 0.0076 + 0.01997 + 0.055 + ... (36) 

The best estimate is formally given by the sum of the first ten terms 

m fix,best = 0.132795 ± 0.000014 (37) 

according to the formal prescription for the the evaluation of precision. The exact 
result m[ = 0.130954, however, does not fall into the tiny interval given by the error 
bars in eq. fl37l) . Therefore the formal criterion of the accuracy is violated in this case: 
the discrepancy rn{ lx,best — m[ = 0.00184 is not controlled by the smallest term of the 
asymtpotic expansion eq. (|3"6"D. Still this discrepancy is small and the actual accuracy 
for the first moment given by the asymptotic expansion eq. (|36|) is 1.3%. This suffices 
for practical purposes. One can easily see the difference between these two observables 
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which reflects the different numerical magnitude of the infra-red contributions. Note 
also that the numerical magnitude of the smallest term of the expansion eq. (^) is 
very sensitive to the value of the third coefficient of the effective /3-function. For the 
decay rate we find r( = 0.1946 and r((l) = 0.1527 to be compared with eq. (p2,23,24|) 
and eqs. (^5|,^6|) . 

^From eq. (|34]) one sees that an infra-red fixed point exists also in fourth order of 
the effective /3-function for any k% < 95.9. For these numerical values the effect of k% on 
the exact moments within the model is rather weak. The pattern of "convergence" for 
the decay rate is mainly determined by the contribution of the zeroth order moment 
(or even by the mixed rh.20 moment) and is very close to the expressions in eqs. (|13|-[35|). 

In this model there are ways of accelerating the convergence with nice results but 
they definitely cannot be justified for use in the general case. Still our conclusion 
about the achievable precision within finite order perturbation theory in fifth order 
remains valid. 



5 Conclusions 

Using the standard estimate of the accuracy of an asymptotic series we have found 
that the theoretical precision in the perturbative description of r-decay is already 
limited by the asymptotic growth of the coefficients in fifth order of perturbation 
theory. This is a scheme invariant statement. The accuracy of perturbative expan- 
sions for a reasonably general set of observables cannot be better than 5% - 10%. 
Taking a stricter attitude we claim that the zeroth order moment is not computable 
within perturbation theory. Any consistent description of r-decay data at fifth order 
of perturbation theory requires exclusion of the zeroth order moment from the list of 
observables (or it should not constitute a dominant contribution). At fifth order of 
perturbation theory and with the present numerical value of the coupling, the first 
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two moments of the spectral density are too different to be simultaneously treated 
by perturbation theory with an accuracy better than 5% - 10%. Therefore one has 
to go beyond finite order perturbation theory to compare these two observables if 
one requires a theoretical accuracy that exceeds present experimental accuracy. This 
implies the use of some procedure of resummation. The resummation procedure is 
not defined uniquely and the result depends on the prescription chosen |16|, |4], [27] . 
Moreover, if one resums the infinite number of perturbation theory terms the con- 
densates have no invariant meaning anymore and their numerical values may change 



24]. Therefore, improving the theoretical accuracy for this system seems to require 
the creation of a new paradigm. 

The extraction of a s from the r-decay rate and its comparison with a s -values 
determined from other experiments does not appear to be the best test of perturbation 
theory for the r system. The crucial test of the applicability of perturbation theory 
for the r system would be the simultaneous calculation of two observables (moments) 
with an appropriate accuracy. If the set of moments includes the zeroth moment 
then the ultimate accuracy of finite order perturbative expansions has been already 
reached. 
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